Quasiparticle scattering by quantum phase slips in one-dimensional superfluids 
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Quantum phase slips (QPS) in narrow superfluid channels generate momentum by unwinding 
the supercurrent. In a uniform Bose gas, this momentum needs to be absorbed by quasiparticles 
(phonons). We show that this requirement results in an additional exponential suppression of the 
QPS rate (compared to the rate of QPS induced by a sharply localized perturbation). In BCS- 
paired fluids, momentum can be transferred to fermionic quasiparticles, and we find an interesting 
interplay between quasiparticle scattering on QPS and on disorder. 
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1. Introduction. Particle production and scattering 
during tunneling are of interest in many applications. 
They can significantly modify the transition amplitude 
and, in addition, destroy coherence between the states 
connected by tunneling. Of particular interest are sit- 
uations when some scattering must occur, as a conse- 
quence of a nonzero momentum generated by the tunnel- 
ing event. 

One of the simplest such cases, and the one on which 
we concentrate in this Letter, is a quantum phase slip in 
a narrow superfluid channel. Advances in containment 
and cooling of atomic Bose gases make it possible to con- 
template a ring (toroidal) geometry — a channel closed in 
the longitudinal (x) direction. This is the case considered 
here, although many of our results are applicable also to 
a cigar-shaped geometry. In addition, since these gases 
are dilute, the correlation ("healing") length ^ can be 
made larger than the transverse size, making the system 
effectively one-dimensional (ID). 

In a ID system, the order parameter ijj can momentar- 
ily vanish at some point, allowing the phase of to un- 
wind. Such events are referred to as phase slips. Here we 
consider the case of very low temperatures when the main 
mechanism for phase slips is quantum, rather than ther- 
mal, fluctuations. These are quantum phase slips (QPS). 
Their high-temperature counterpart, thermally activated 
phase slips (TAPS), was considered in the framework of 
the nucleation theory of Refs. 0,0,0 in Ref- B- 

The purpose of this Letter is to present a calculation 
of how the QPS rate is affected by interactions of QPS 
with quasiparticles (phonons). First, we show that unless 
there is a highly localized (or very strong) external per- 
turbation, QPS without momentum transfer to phonons 
are strongly suppressed by destructive interference. This 
makes such QPS highly unlikely and makes interaction 
with phonons a principal effect. Next, we compute the 
rate of QPS assisted by phonon scattering and find that 
at low temperature T it is given by 7?. '--^ exp{—VsP/2T), 
where P — 2'nn is the momentum generated by unwind- 
ing the phase; n is the linear number density of the gas, 
Vs is the speed of Bogoliubov's phonons. We also es- 
timate the preexponent. The physical interpretation of 
this result is that the destructive interference is removed 
when momentum P is absorbed by the phonon subsys- 



tem, but that requires energy of at least VsP/'i (to allow 
for tunneling between phonon states with total momenta 
-P/2 and P/2). 

The above estimate for the rate is the main result of 
the present paper. It applies to a uniform weakly-coupled 
Bose gas at temperatures T <^ gn and represents a strong 
suppression, compared to the case when QPS is induced 
by a highly localized perturbation, and momentum con- 
servation plays no role. We briefly discuss means to ex- 
perimentally detect phonon scattering on QPS in atomic 
superfluids. Finally, we compare this case to the case of 
BCS paired superconductors. 

2. QPS without momentum transfer to phonons. By 
unwinding the supercurrent, each QPS releases momen- 
tum ±P. (In the ring-shaped geometry, a better con- 
served quantity would be the angular momentum, but 
this distinction will not be important in what follows.) 
This corresponds to an imaginary contribution to the Eu- 
clidean action of a QPS: 

^'S'qps = —2Tri / n{x)dx « —2TTinxQ , (1) 
Jo 

where xq is the QPS location, and n is an average density. 
Then, the tunneling amplitude is 



A— I dxoB{xo)e 



(2) 



If the scale of variation of B{xo) (due to possible vari- 
ations of density, etc.) is some large £, the ampli- 
tude 1(21 is suppressed by exp(— const. n£). For £ ^ 
this suppression is much stronger than the semiclassical 
exp(— const. Vg/g). It can be seen as a consequence of the 
destructive interference among QPS with different values 

of Xq. 

One way to avoid this suppression is to consider a lo- 
calized perturbation with scale £ -^i £■ Calculations for 
this case have been done in Ref. [5|- Another possi- 
bility, which we consider here, is that the momentum 
released by a QPS is absorbed by phonons. We expect 
that if phonons change their total momentum by 27rn, 
the recoil on the QPS will be zero, and the suppression 
by exp(— const. n£) will be lifted. For sufficiently uni- 
form gases, such phonon-assisted tunneling may be the 
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dominant QPS mechanism. Furthermore, the two mech- 
anisms can be distinguished experimentally by measuring 
the momentum transfer (see below). 

3. Phonon production at T = 0. In this case, the 
phonon subsystem is initially in the ground state, so to 
transfer momentum to phonons we need to produce them 
in the final state. To set up come notation, we briefly 
discuss this case and show why, if the only source of en- 
ergy available for phonon production is the energy due to 
unwinding of the supercurrent, such processes are kine- 
matically forbidden. The formalism sketched here may 
be useful when there are other sources of energy, as for 
instance in various "cosmology in the lab" proposals @, 
where the properties of the superfluid are changed in time 
to mimic certain cosmological spacetimes. 

We consider a weakly-coupled ID Bose gas. In ID, the 
effective coupling constant g has dimension of velocity, 
and the dimensionless measure of the coupling strength 
is 



— = 0.008 



VlOnm/ 



/2 /IQi 



where Vg in the phonon speed, a is the scattering length, 
n'^'^'> is the 3D number density, and R is the transverse 
radius of the channel. As we will see, the real part of the 
instanton action is of order irvg/g (times a logarithm). 

The crossover between the weakly and strongly coupled 
regimes can be located using the results of Haldane : 
it occurs at g/vs ~ tt. In what follows we assume a 
somewhat weaker coupling g/vs ~ 1, which according to 
(O can be achieved for i? ~ 0.1 /im. Note that such 
small transverse radii may be necessary in any case, just 
to drive the system into the ID regime. Indeed, the ID 
regime requires i? <C while the typical values of ^ are 
a few tenths of micron. 

So, if it were not for the destructive interference, an 
instanton action of order S'qps ~ T^Vs/g would allow for 
a relatively large QPS rate in narrow channels. 

Since the system is in the weakly coupled regime, we 
can obtain a more precise estimate for S'qps by using the 
Gross-Pitaevsky (GP) description. This is given by the 
Euclidean Lagrangian 

Le = i^^dryj + ^IS.V'P + f IV'I' - MIV-P . (4) 

where t = it is the Euchdean time. Outside the QPS 
core, fluctuations of density are small, and we can lin- 
earize the system by writing ip = \/n + Sn exp{i9) and 
expanding in small Sn. Integrating out Sn and assuming 
the long- wavelength limit, we obtain 



indrO + — {drOf 
2.9 



2M^ ' 



(5) 



From this, we read off the sound velocity Vs = ygnjM . 
In what follows, we assume that n is x-independent and 
choose our unit of length so that = 1. 



Since (0) is a long- wavelength limit, it implies a cer- 
tain choice of the reference frame: momenta of produced 
phonons should be small in comparison with 1/^. Thus, 
the Galilean invariance of the original theory Q now 
holds only approximately and only for boost speeds much 
smaller than (7\f^)^^ ~ Vg- This approximate symmetry, 
however, would be completely sufficient to restore the 
topological (first) term in eq. ©. 

Instantons corresponding to QPS are vortices of the 
theory in the (a;,r) plane. Using the limit one 
finds that away from the core the phase 9 of the instanton 
solution is 

9i{x - xo, r - To) = arg[a; - xq + «(t - tq)] , (6) 

where (xq, tq) is the location of the QPS. 

Consider an instanton-antiinstanton pair in the pres- 
ence of a persistent supercurrent. The corresponding 
classical configuration is 

27r 

9c{x,t) = 9i{x - xq.t - To) - 9i{x,t) + -j^Nx (7) 

(here we have placed the antiinstanton at the origin). 
The last term corresponds to a supercurrent with winding 
number N over the longitudinal length L of the system. 
Substituting ||7J) into Q and integrating, we obtain the 
Euclidean action of the pair to logarithmic accuracy 



Spair = ~2TrinxQ — Etq H In - + O ( - 

9 ^ \9 



(8) 



where d ^ {xl + t^Y^^, and ^ = {AgMn)-^/^ is the 
healing length (the size of the vortex core). Eq. © 
applies when the separation d is much smaller than the 
longitudinal size L of the channel. Contribution from the 
QPS core is subleading and included in the 0{l/g) term. 

The first term in |(SJ) is the imaginary part anticipated 
in QJ. The second term, with 



E ^ E 



N 



(9) 



comes from the cross term between the instanton and 
the supercurrent contributions to dx9. Note that E^ is 
precisely the energy released by a single phase slip. 

It is easy to see, however, that this energy alone is not 
sufficient to satisfy the kinematic condition for produc- 
tion of phonons with total momentum P = 2Tm. Indeed, 
the minimal (critical) energy required for that is 



Ecr ~ 2'Kn = — 
.9^ 



(10) 



Landau's criterion for superfluidity is N/L < I/Att^, so 
in the superfluid state En < Ecr [En = Ecr precisely at 
the critical current). We conclude that at zero tempera- 
ture QPS are effectively impossible, unless one considers 
short-scale perturbations or introduces other sources of 
energy. This conclusion is consistent with Galilean in- 
variance (GI) of the T = state. One way to break GI 
is to include thermal fluctuations. 
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4- Phonon scattering at T ^ 0. In this case, there 
are preexisting phonons in the system, and we consider 
a QPS assisted by scattering of phonons, so that the mo- 
mentum released by the QPS is absorbed by the phonon 
subsystem. We take the hmit when the energy is 
neghgible compared to the energy of the requisite phonon 
state. At T ^ 0, the system can be regarded as compact- 
ified on a cyhnder of circumference (3 = 1/T. Equiva- 
lently, instead of a single lA pair we should consider a 
periodic chain of such pairs: 

0c = arg[l - e-'^r-ro)-{a:-Xo)^ _ ^j.g(^ _ g-^r-x^ _ (^^^ 

Here, instantons are shifted relative to antiinstantons by 
amount xq in the x-direction, and by tq in the r-direction. 
These are "soft" collective coordinates, which will need 
to be integrated over. All lengths in Hll() are in units of 
/3/27r, and we assume that the total longitudinal length 
L is effectively infinite. 

The action per period, to logarithmic accuracy, equals 

2tT 1 TT 

Schain ~ ^"iP^Q-l In 1 ln(cOSh — COS Tg ) . (12) 

The last term is the instanton-antiinstanton interaction. 
If it were not for the momentum P, the dominant config- 
uration would have xq — and coincide with the periodic 
instanton which describes tunneling between quasi- 
particle states with zero total momenta. At nonzero P, 
configurations with xq ^ become important. 

The integrals over xq and = —itq can be done by 
steepest descent. In the limit /? 3> C (that is, T ^ gn), 
the saddle point is at tq — t: and xq ~ iir. The 
saddle-point action is dominated by the first term in (|12|l : 
•Schain ~ P/2T (where we have returned to the physical 
units of length). Corrections to this result are controlled 
by the small parameter T^. 

The above action determines the exponential factor in 
the QPS rate. The preexponent is given by the determi- 
nant of field fluctuations near the saddle point. As usual, 
the main effect is due to the collective coordinates. Two 
zero modes, corresponding to translations in space and 
time, each contribute a factor of order l/^y^, up to a 
power of ln(/3/^). The saddle-point integrals over xq and 
To result each in a factor of order (A:*^)~^/^, where fc* is 
the typical momentum of phonons scattered by the QPS. 
Other modes result in a factor of order unity. Assembling 
all the factors together, we obtain, for the QPS rate per 
length L, in the hmit (3 ^ ^, 



g^^k^ 

(where we have restored Vs). A detailed study of the 
initial and final phonon states (to be presented else- 
where) gives kl ~ (^^/?)~^ ln(/3/^), in the same limit. 
So, fc, ^ which confirms that the longwave limit 
is applicable to the present problem, with corrections 
controlled by fc,^. 



The exponential in (|13|l can be interpreted by noting 
that at a finite temperature the quasiparticle system can 
tunnel from a state with momentum — P/2 to a state 
with momentum P/2, thus absorbing the full momentum 
P but requiring energy of only P/2. The main physical 
effect reflected in (|13|) is that transfer of a large momen- 
tum requires tunneling from a fairly high-energy state, 
resulting in an additional suppression. If no momentum 
transfer were necessary, the leading term in the expo- 
nent would be {2TT/g) In(r^) ^ (see also Ref. so the 
suppression would be much weaker. 

Phonon scattering by QPS in atomic superfluids can 
in principle be detected experimentally by measuring the 
momentum distribution of atoms, using, say, the stan- 
dard momentum imaging. In the ring-shaped geometry, 
this can be done by opening the trap in one place and 
looking at the atomic cloud emitted from there. For ex- 
ample, if we start with the state without supercurrent, a 
single QPS will create a unit of supercurrent and a com- 
pensating normal flow in the opposite direction. Since, 
the "normal" atoms have comparatively large momenta, 
the cloud will be highly asymmetric. 

5. Quasiparticle production in superconductors. The 
case of BCS superconductors is especially interesting be- 
cause nonzero resistance has been observed in ultra-thin 
wires at low temperatur es |M llC| , and it has been inter- 
preted in terms of QPS jllf . 

Curiously, the "hydrodynamic" part of the effective 
Lagrangian for superconductors is given by the same ex- 
pression Q as for atomic superfluids, except that ip is 
now the field of Cooper pairs, which describes charge 
fluctuations in the wire (cf. Ref. Accordingly, the 

coupling constant is g = 4e^/C, where C is the capaci- 
tance per unit length. The reason for this similarity is 
that in ID screening is weak, and the plasmon mode is 
gapless (the Mooij-Schon mode 0|). From 10}, we ob- 
tain the speed of plasmons as cq = (e^ns/m,C)^/^, where 
Ug = A{'ip'^'4')m^/M, and to* is the electron mass. This is 
indeed the speed of the Mooij-Schon mode. 

For BCS superconductors, there is an additional La- 
grangian describing the interaction of if) with fermionic 
quasiparticles. Indeed, production of these quasiparticles 
is the least energy-consuming way of momentum produc- 
tion, since each quasiparticle carries a large momentum 
kp but costs only a relatively small energy equal to the 
superconducting gap A. This mechanism is unavailable 
to atomic superfluids. 

On general grounds, one expects that impurities, which 
break both the Galilean invariance and momentum con- 
servation, may significantly affect the physics of QPS [l3 |. 
Here we present an explicit example of such interplay, 
using the simple, if somewhat artificial, single-channel 
limit, in which electrons are confined to motion in the x 
direction only, and their mean-free path A is taken to be 
much larger than the QPS core size. 

The Fermi surface now consists of two points, at fc = 
ifc^, so there are right-movers, denoted by subscript (-I-), 
and left-movers, with subscript (— ). Each variety can 
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have up or down spin, corresponding to a =t,|. We 
define new fermionic operators a„± related to the original 
electron operators Ca± by 

c,±(x,i) =e^'''""a.±(a:,i) . (14) 

The new operators can be conveniently assembled into 
two Dirac spinors: — (a|+,a|-), — (a|_|^,a|_). 

In the case without disorder, A oo, the Lagrangian 
of fermions can be written in the relativistic form: 

Lf^iY. *-7^5^*. - 5V^T*i - 5V*^i*T ■ (15) 

Here = 0, 1, the (1+1) 7-matrices are 7" — (Ti, 7^ — 
—i(T2, and — 5*^70 • In this section, we set vp = 1- 

For any ip with a nontrivial winding at infinity, the 
Euclidean equation of motion for 'J has a nontrivial nor- 
malizable solution (zero mode) . The zero mode is known 
explicitly for -0 of the form ijj{r, -d) — f{r) exp(i?9) (in po- 
lar coordinates) ^15] , and its existence in the general case 
is guaranteed by an index theorem . The zero mode 
causes the amplitude for tunneling without fermions to 
vanish identically. As has been explained by 't Hooft 
\it\ in the context of gauge theories in four dimensions, 
this means that only certain anomalous Green functions, 
involving fermions, are nonzero. In our case, the rele- 
vant Green function is (a|_(x)a|+(2/)). It corresponds 
to creation of a quasiparticle-quasihole pair at the QPS 
core. The total momentum produced in this process is 
2kp. The phenomenon is a ID analog of "momentogen- 
esis" 01 by vortex motion in a 2D array of Josephson 
junctions. 

Note that 2fci? is precisely the momentum necessary 
to cancel the momentum —P released by unwinding of 



the current. Indeed, P = 2T:{tp'^^) = nn, where n is the 
electron density. Since n — 2kp/TT, we have P — 2kF- If 
the spatial dependence of the zero mode for a|+ is x{x — 
xq), then for the original operator c^^ it is eKp[ikp{x — 
xq)]x{x — xo). In this way, we find 

(c|_(x)cT+(y))^ 

J dxodroe-^'e'^-(^+«)x(x-xo)x(y-xo) , (16) 

where x = (x, r), and S' is the real part of the instanton 
action. 

At r = and in the absence of disorder, production of 
a quasiparticle-quasihole pair via the matrix element 
is kinematically forbidden in the same way as phonon 
production was in the case of an atomic superfluid (Lan- 
dau's criterion now being that the velocity of the super- 
flow must be smaller than the depairing velocity A/kp). 
However, scattering on an impurity potential W{x) in- 
troduces an additional Lagrangian 

ALp = W{x)Y,ci-Ca++ii-c. . (17) 

(7 

Computing to the first order in W and absorbing the 
fermion operators into the anomalous correlator (|16|l . we 
find that in the presence of disorder the amplitude of 
QPS without quasiparticle production becomes nonzero. 
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